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Kelly-Borges et al. (1991) and Kelly-Borges & Pomponi (1994) utilised partial 18S rRNA 
gene sequences to resolve relationships within hadromerid and lithistid sponges (Porifera: 
Demospongiae). While their results clarified several specific systematic problems, their 
conclusions were hampered by low levels of sequence variation. This study sought primarily 
to evaluate the resolution potential between regions of the 18S rRNA gene used in previous 
studies on sponges, and 28S rRNA genes used in more recent work. Six lithistid sponge taxa 
were chosen to represent a gradient of taxonomic relationships, ranging through genus, 
family, order and class. Approximately 1,300bp of the 18S rRNA gene and a 700bp region at 
the 5’ end of the 28S rRNA gene were compared with the data of Kelly-Borges & Pomponi 
(1994). We found that the 700bp region of the 28S rRNA gene presented the greatest 
potential for resolution of this group of Porifera at the genus and family level, and that the 
resultant molecular phylogeny was congruent with morphological hypotheses for the group. 
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For the past twenty years, organismal phylo- 
genies have been inferred from the primary 
sequence ofa portion of their genomes. The small 
subunit ribosomal RNA (SSU rRNA) or 
18SrRNA gene has dominated the field as 
molecule of choice. The analysis of this molecule 
has been singularly instrumental in elueidating 
the phylogenetic relationships and natural history 
of almost all known prokaryotic species where 
attempts using other methods have failed (see 
Woese, 1987). Following the success of 
microbiologists adopting this approach, the 
systematics of eukaryotic taxa has been addressed 
by sequencing the 18S rRNA gene (Sogin et al., 
1986), This gene is probably still the most 
frequently used for this purpose. 


It is often desirable to use more than a single 
gene region for the reconstruction of a phylogeny, 
to supply additional, potentially eorroborative 
phylogenetie hypotheses. For example, the 5’ 
region of the 28S rRNA gene has also been used 
with great effeetiveness (Baroin et al. 1988; 
Chombard et al., 1998), as have Elongation Factor 
genes (Iwabe et al., 1989; Rivera & Lake, 1992), 
ATPascs (Iwabe et al., 1989) and DNA- 
dependent RNA polymerases (Puhler ct al. 1989), 
among others. There are a number of additional 


prerequisites for choosing a gene for the purposes 
of reconstructing a phylogeny. It is desirable that 
a eonstaney of funetion (funetional orthology) is 
maintained throughout the evolution of the taxa 
of interest. Problems associated with long 
branches may be observed on trees where some 
genes have experieneed a relaxation in sclective 
pressures. The possibility of mistakenly isolating 
a paralogous homologue must be minimal. A 
suitable genc must show signs of having enough 
variability to discriminate between taxa at the 
desired taxonomic Icvel. It must also be 
conserved enough to permit robust alignments 
and comparisons across the deepest divisions. 


Sponge phylogenies, in particular, have given 
rise to a number of contentious arguments, most 
of which result from a laek of suitably variable 
morphologieal characters which distinguish 
sponges at the species level and higher (Van Soest, 
1987; Hadju et al., 1993; Kelly-Borges & 
Bergquist, 1997; Sandford & Kelly-Borges, 
1997). The primary diagnostic morhological 
eharacters that differentiate sponge genera and 
species are spicule morphology and their 
arrangement within the sponge body. Characters 
often influeneed by environmental factors, such 


as texture, surface features and colouration are 
less reliable as they are frequently plastic. 


Although the construction of poriferan 
phylogenetic hypotheses using molecular 
sequence data is still at a preliminary stage, with 
very few studies completed, it is likely that many 
future sponge systematic projects will incorp- 
orate a molecular moiety. For this reason, it is 
necessary to establish the taxonomic levels at 
which certain gene regions will be appropriate. 
For example, some stretches of DNA might be 
informative about phylogenetic relationships at 
the genus and spccies level, but they may be 
unsuitable tor studies of ordinal relationships and 
so on. All studies so tar have used different genes 
and gene regions for phylogenetic purposes 
(Kelly-Borges et al., 1991, 1994; West & Powers, 
1993), preventing any useful links between thesc 
data towards the construction of larger phylo- 
genies. The advantage that may be gained in 
future by using the samc gene region in all studies 
is therefore obvious. 


The primary goals of our study were to evaluate 
the resolution potential between various regions 
of sponge 18S and 28S rRNA genes, in order to 
determine which genes would efliciently resolve 
phylogenies at several taxonomic levels. Because 
of past difficulties in resolution using [8S rRNA 
(Kelly-Borges et al, 1991; Kelly-Borges & 
Pomponi, 1994), we took a positive approach in 
our more recent research to determine which gene 
would successfully provide resolution within a 
group of lithistid sponges, and thus, potentially 
within other taxonomic groups. 


In this study we evaluate the relative utility of 
four alignments with diffcrent gene origin, 
sequence length, and method of analysis. To do 
this, we extended the 18S rRNA gene data of 
Kelly-Borges & Pomponi (1994) for six specics 
up to approximately 1,300bp, and in addition, 
using the same taxa, we have sequenced approx- 
imately 700bp of the 5’ end of the 28S rRNA 
gene. Taxa were selected to encompass a range of 
taxonomic levels including genus, family, order, 
and class. The criteria by which the sequence data 
sets were evaluated for their potential utility 
included counting the number of variable sites 
and the number of parsimony-informative sites, 
using maximum likelihood in order to estimate 
the proportion of constant sites that might be 
invariable, and also conducting an objective analysis 
of the resulting tree topologies. The resulting 
topologies were compared with a hypothetical 
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reconstruction based upon morphological 
characters. 


For the purpose of this exercise four lithistid 
sponges (Class Demospongiae) were selected 
from a much larger study on sponge phylogeny, 
and two hexactinellid sponges (Class Hexact- 
inellida) were chosen as an outgroup. Lithistid 
sponges represent relict forms of an ancestral 
fauna from which, it is thought, most demo- 
sponges have evolved. These sponges are 
characterised by the possession of a rigid siliceous 
skeleton made up of irregular branching desma 
spicules, the ends of which interlock (zygose) 
with neighbouring spicules (see Kelly-Borges & 
Pomponi, 1994, Figs 1,2). In some cases there are 
additional spicules present, providing clues as to 
their affinities with non desma-bearing sponges, 
and the polyphyly of at least some of these genera 
has been recently confirmed by Kelly-Borges & 
Pomponi (1994). Although very difficult to 
differentiate morphologically, the most reliable 
diagnostic characters that can be used are the 
morphology, ornamentation, and pattern of 
zygosis of the desma spicules, and the morpho- 
logy of the additional spicules ifthey arc present. 


MATERIALS AND METHODS 


TAXA SELECTION. Sponge species were 
selected from a broader study of sponge phylo- 
geny, specifically for the purposc of examining 
capability of taxonomic resolution of poriferan 
sequence data (Table 1). All sponges were 
collected using Harbor Branch Oceanographic 
Institution’s ‘Johnson-Sea-Link IP manned 
submersible, except Theonella spp. which were 
collected using SCUBA. Samples were identified 
through histological examination of skeletal 
structures, the procedures for which are detailed 
in Kelly-Borges et al. (1994). Voucher specimens 
have been deposited in the collections of either 
The Natural History Museum, London (BMNH), 
or the Harbor Branch Oceanographic Museum, 
Florida (HBOM). Registration numbers are given 
in Table 1. Hexactinellid sponges Margaritella 
coeloptychioides and Sympagella nux (Table 1) 
were selected to provide outgroup sequences for 
the molecular phylogeny reconstruction. 


MORPHOLOGICAL PHYLOGENY RECON- 
STRUCTION. A phylogenetic analysis of 
morphological characters (Table 2) was carried 
out to examine relationships within and between 
the theonellid and corallistid taxa for the purpose 
of comparison with trees gained from recons- 
truction of molecular data. The analysis used the 
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TABLE 1. Collection data and taxonomic position for sponge 


taxa sequenced in this study. 


Taxon 


Museum 
Registration 


Locality 


Class Hexactinellida _ 


Subclass Hexasterophora 


Order Hexactinosida 


Family Euretidae 


Margaritella coeloptychioides 
| Schmidt, 1870 


HBOM 003:00925 


Turks and Caicos 


Order Lyssacinosida 


Family Caulophacidae 


Svmpagella nux 
Schmidt, 1870 


HBOM 003:00929 


Turks and Caicos 


Class Demospongiae 


Subclass Tetractinomorpha 


| Order Astrophorida 


Family Theonellidae 


| Theonella sp. 1 


BMNH1998.3.4.1 | 


Belau, Micronesia 


Theonella sp. 2 


BMNH1998.3.4.2 


Belau, Micronesia 


Discodermia sp. 


BMNH1998.3.4.3 


Bahamas, Caribbean 


Family Corallistidae 


Corallistes tvpus 


Schmidt, 1870 


BMNH1998.3.4.4 


TABLE 2. Morphological characters (A) and their 
character-states (B) for taxa in this study (see Kelly-Borges & 
Pomponi, 1994, for an explanation of characters). Characters 
indicated as * are absent from the outgroup in that they do not 


possess desmas. 


Branch and Bound search option of PAUP 
3.1.1, and data were unordered and 
unweighted. In order to obtain a dirccted 
analysis, members of two non-lithistid 
astrophorid families, Geodia (Family 
Geodiidae) and Stelletta (Family Ancor- 
inidae) were chosen as outgroups. The 
major skeletal charactcrs that separate 
these lithistids from their outgroups are 
the possession of desmas, unique 
ornamented dichotriaenes, and certain 
types of microscleres (see Kelly-Borges & 
Pomponi, 1994). 


MOLECULAR EXPERIMENTAL 
PROCEDURES. Sample collection, 
prescrvation and DNA extraction have 
been previously described (Kelly-Borges 
& Pomponi, 1994). PCR primers and 
sequencing oligonucleotides are listed in 
Table 3 for both genes. PCR reactions 
were carried out ina 50l reaction volume 
which contained a one-tenth volume of 
10x PCR buffer (500mM KCI, 100mM 
Tris-HCl, 1% Triton X-100), dNTPs to a 
final concentration of 200uM, primers at a 
concentration of 200uM and 2.5mM 
MgCl. The PCR protocol began with an 
initial denaturation at 94°C for 5mins, 
followed by 35 cycles of denaturation at 
94°C for | min, annealing at 55°C for Imin 
and extension at 72°C for 1min. This PCR 


A, Character x i 5 
number Character Character state regime was used for both genes. Following 
i Desmas: a, present; b, absent cycling, the success of the amplification 
= Desma a, tetraclonal; b, dicranoclonal; was determined by electrophoresing on an 
= development. c, absent* ethidium bromide-stained agarose gel and 
a, articulated at ends of zygomes visualised by short-wave UV illumination. 
Zygosis (Fig. 2A,C); For each taxon, a total of 10 PCR reactions 
3 : b, articulated along zygomes ` 
architecture. (Fig. 2E); were carried out and the products were 
c, absent* pooled. This was an effort to reduce the 
tinal a, oxea; ; potential for an amplification-induced error 
i; megascleres. b: a ea in the sequences. The pooled amplification 
i TE products were electophoresed on a single 
a, short-shafted discotriaenes l d the band ised 
(Figs 1C, 2B); agarose gel and the band was excise 
i : THaene ee phyllotriaenes using a clean scalpel blade. The DNA 
f ig. 2D); 
megascleres. c, long-shafted ornamented dicho- and 
| trichotriaenes (Fig. 1A,B); 3 
D d, long-shafted ortho- and dichotriaenes | || B- Taxon TER We A S6 i S 
ji 6 Euasters. a, present; b, absent | Discodermia aļ| ajaja | a| bl[bljala 
Theonella sp. | al b]b/| bib b 
T eral a, present; b, absent eee aja z 
ll (Fig. 2F). Theonella sp.2 |a | a |ajb|bjb|b]ajļb 
Small acanthose li : z b b blb 
8 microrhabds a, present; b, absent Sael tsres SPs 2 a = a 
(Fig. 2B). | | Geodia Ble lela lala cb ble 
(outgroup) 
Large acanthose oe L 
9 microrhabds a, present; b, absent telletta E 
(Fig. 2B). (outgroup) bye jecjajdjajbjbjb 
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fragment was extracted from the agarose usmy 
the Quiaex If (Quaigen Lid, UK.) PCR 
purification kit. DNA sequencing reactions were 
varied out using the Amplitag FS sequencing kit 
{Applied BioSystems, Inc.). The sequencing 
protocol for both genes was carried out as per the 
manufacturer's instructions. Sequencing was 
carried out on an Applied Biosystems 373 
automated sequencing apparatus and the data was 
analysed using the Sequence Navigator software 
(Applied Biosystems Inc.). We estimate that 98% 
of gene regions were covered by more than one 
sequencing read, with approximately 25% being 
covered by three or more sequencing overlaps. 
‘The sequences haye been deposited in the EMBL 
sequence repository under the accession numbers 
AJ224646-AJ224651 (demosponge 185 rRNA). 
AJ224123-AJ224] 24 (hexactinellid 185 rRNA), 
AS00591 -AJ005918 (288 rRNA), 


MOLECULAR PHYLOGENY RECON- 
STRUCTION. Nucleotide positions whose 
identities were not possible to establish unamb- 
iguously, Were coded according to the 
International Union of Pure and Applied 
Chemistry (TUPAC) nomenclature. The sequences 
were aligned using the Genetic Data Environ- 
ment (Smith, 1993). In the majority of cases, the 
positional orthology of the nucleotides was 
relatively casy to establish. A conservative 
approach to the alignment was taken, with those 
positions whose homology was net possible to 
establish with absolute certainty, being excluded 
from subsequent analyses, Attempts in the most 
difficult cases, to relate the Sequences to each 
other, on the basis of RNA secondary structure 
proved recalcitrant, In the absence of conclusive 
grounds for establishing homology, the sites in 
question were excluded, Phylogenetic hypothesis 
construction and sequence statistics were 
evaluated using PAUP*4.0d54 (Swofford, 1993). 
Transition-transversion ratios were calculated 
from each dataset by first constructing a 
neighbor-joining iree from LogDet distances 
(Lockhart etal., 1994), This tree was considered a 
working hypothesis of relationships, The 
estimate of transition-transversion ratio might he 
influenced to some extent by tree topology, but 
this influence was not thought to be significant. 
Using maximum likelihood criteria, the 
transition-transversion ratio was chosen that 
yielded the highest likelihood value. The gamma 
shape parameter fùr each dataset was calculated 
using the optimised transition-transversion ratio, 
again using maximum likelihood as the 
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optimisation criterron. The samma shape param- 
eter that yielded the highest likelihood was 
chosen. 


RESULTS 


MORPHOLOGICAL PHYLOGENY RECON- 
STRUCTION, A single minimum length tree 
was obtained of length 14 anda very high consist- 
ency index (Cl) of 1.0. The phylogenetic tree 
hypothesises that species of Theonella are more 
recently derived (han Diseodermia, and that these 
two genera form a clade more recently derived 
than Curallistes. This topology is identical tu one 
ef the reconstructions derived from sequence 
data (Fig. 3A), and supports the current classif- 
ication that recognises the differentiation of 
Discedermia and Theonella ftom Corallistes in 
the Family Theonellidae and Family Corallist- 
idae, respectively (see Kelly-Borges & Pomponi, 
1994), 

Morphological characters 2,3,7, 8 and 9 (Table 
2) differentiate Diseoderniia and Theonella from 
Crrallistes, and the states of characters 4 and $ 
differentiate Theonella from Discodermia, and 
both from Corallistes. In Discodermia and 
Theonella the desmas are tetracrepid (character 
2a) with four clones (zygomes) which clasp 
(zygose) at their very ends (character 3a; Fig. 
2A4,C). Corullistes, on the other hand, has dicrane- 
clonal desmas (character 2h) in which zygosis 
spreads. along the zygomes (character 3b; Figs 
1A, 2B) 

The desma skeleton of these genera is 
supplemented with monaxonal megascleres - oxea 
in Discodermia (character 4a), long and curved 
with blunt hammer-like ends (tylostrongyles) in 
Theonella (character 4b; Fig. 2C), and whispy 
roughened oxea-like spicules in Corallistes 
(character 4c, Fig. LA). Triaenose megascleres 
are found at the surface of the sponge with their 
head rays parallel with the surface and the rhabd 
perpendicular to the surface (Fig. 1A,C). These 
are discotriaenes in Discodermia (character Sa; 
Fig. IC, 2B), phyllotriaenes in Theonella 
(character 5b; Fig. 2D), and ornamented 
(character Sc: Fig. LAB) or plain (Fig. 2F) 
dichotriaenes in Corallistes. Microseleres are 
scattered throughout the sponge body and otten 
form a thick surface erust, In Discodermia here 
are lwo size categories of roughened 
microrhabds (characters 8a, 9a; Fig. 2B), 
whereas in Theonella there is onty one (character 
Ba; Fig. IC). Corallistes lacks microrhabds but 
possess sireptasters (Fig. 2F). 
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TABLE 3. Oligonucleotide names and their 
corresponding sequences. The first and second 
oligonucleotides are designed to amplify a large 
portion of the 18S rRNA gene, whilst the third and 
fourth oligonucleotide sequences are designed to 
amplify an approximately 700bp stretch of the 28S 


tRNA gene. o 

| a a Sequence | 
TGG TAC GGT AGT GGC CTA CCA TGG 

[ACG GGC GGT GTG TAC AAA GGG CAG | 

GAC CCG TCT TGA AAC ACG A | 

ACA CAC TCC TTA GCG GA | 


|__ 185120 
| 188421 
i RD3A 


RD5B2 


MOLECULAR PHYLOGENY RECON- 
STRUCTION. In total, six taxa had a portion of 
their 18S genes sequenced. The final alignment 
was in excess of 1,300bp in length. The positions 
whose alignment could not be determined unamb- 
iguously were removed. When all taxa were 
considered, the final alignment was 1,135 
positions in length, whereas with the exclusion of 
the hexactinellid sequences, the number of align- 
able positions increased to 1,269. For the 28S 
rRNA gene dataset, almost 700bp were 
scquenced for each taxon. The final alignment 
lengths were 505 positions for the eight taxon 
dataset and 583 positions for the ingroup taxa 
alone. 


A total of eight alignments were analysed for 
their information content. Statistics that were 
evaluated included the number and percentage of 
variable sites, the number and percentage of 
parsimony-informative sites, the estimated 
gamma shape parameter for rate variation across 
sites and the transition-transversion ratio. The 
results of these analyses are given in Table 4. The 
top half of the table refers to the alignments that 
were used when the two hexactinellid outgroup 
sequences were included, whereas the bottom 
half of the table refers to the ingroup only (in this 
instance lithistid demosponges only). 


In most cases, the exclusion of the 
hexactinellid outgroup sequences facilitated the 
use of a longer gene region. This was due to the 
difficulty of aligning some hexactinellid regions 
with their equivalent location in the demosponge 
genes. On the other hand, removal of the 
hexactinellid sequences had the effect of 
reducing the numbers of variable and informative 
sites, with a consequent increase in the number of 
constant sites. The combined dataset always 
contained the highest number of variable and 
informative sites as would be expected, but in 
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FIG. 1. Skeletal architecture of lithistid demosponges. 
Transverse sections have been taken through the 
sponge surface, cellular material removed using 
hydrogen peroxide, and viewed by SEM. A, 
Corallistes nolitangere. Transverse section through 
surface of sponge showing dicranoclonal desma 
reticulation (D). Long-shafted ornamented dicho- 
triaenes (T) emerge from the desma reticulation with 
the shaft perpendicular to the surface and the 
cladomes (head) parallel with the surface. Oxeote 
spicules (O) resembling fine hairs can be seen in 
residue cellular material. Streptaster microscleres 
pack the surface (M). Scale=238um. B, Corallistes 
nolitangere. View of the sponge surface showing the 
ornamented heads of the dichotriaenes, with the desma 
reticulation visible beneath. Scale=400um. C, 
Discodermia sp. Transverse section through surface 
of sponge showing tetracrepid desma reticulation (D) 
and above this a dense crust of two sizes of acanthose 
microrhabd microscleres (M). Short-shafted disco- 
triaenes line the surface with discs overlapping. 
Scale=11 lum. 


neither instance did it contain the greatest perc- 
entage. One of the striking features of both kinds 
of analysis is the performance of the 28S rRNA 
gene dataset. This region had the highest 
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HIG. 2. Desma, triacne and microsclere morphology in Discodermia, Theonella and Corallistes. A-B, 
Discodermia sp. A, Tetracrepid desmas with tuberculate zygoses (Z) at the ends of the zygomes. Scale=125um; 
B, Short-shatted discotriaene (D), large and small acanthose microrhabd microscleres (M), and large regular 
oxca(O). Scale=20um. C-D, Thenonella sp. C, Tetracrepid desmas with zy goses (Z) at the ends of the zygomes. 
Acanthose microrhabd microscleres are of one size (M) and strongyles have tylote ends (S). Scale=30um. D, 
Short-shafted phyllotriaene. Scale=50um. E-F, Corallistes nolitangere. E, Dicranoclonal desmas with zygoses 
(Z) along the zygomes. Scale=100um; F, Long-shafted dichotriaene (D), streptaster microscleres (S). Scale=50um. 


RESOLUTION POTENTIAL OF 18S AND 28S rRNA SPONGE GENES 349 


TABLE 4. Results from analyses of cight alignments. The top half 
of the table refers to the alignments that were used when the two 
hexactinellid outgroup sequences were included, whereas the 
bottom half of the table refers to the ingroup only (i.e. lithistid 
demosponges only), where Sympagella nux and Margaritella 
coeloptychioides were excluded. The first column contains the 
gene region that was used for each particular analysis. The words 
‘all taxa’ in parentheses indicates that all six taxa were used in that 
particular analysis. 18S (short) refers to the data of Kelly-Borges 


eight taxon datasets indicates that the 
addition of more sequences caused more 
variation at sites that were already 
variable, whilst conserved sites 
remained so even with the addition of 
more distant taxa. 


Of the twelve possible substitution 


types, there are twice as many possible 
transversion mutations as transition 


& Pomponi u 994). 


Length Informative 


5 maa 
Gamma Ti/Tv | 


mutations. For these datasets, the 


Gene region Variable 2 Ben ae eas 
| 18S (long) (all taxa) | 1138 | 446 (39%) | 163(14%) | 0.31 0.98 maximum likelihood transition-trans- 
: > — 5 l | version ratio is slightly above 1.0 in 
| 188 (short (all taxa) | _473_| 137 (29%) | 75 (16%) | 0.24 | 1.20 most cases (except for the 18S rRNA 
_ 288 (alltaxa) | 505 _| 142.(28%) | 84(17%) | 0.24 | 179 | dataset with the eight taxon 
Total (all taxa) 1643 | 588 (36%) | 247 (15%) 0.29 | L18 alignment). This is indicative of a bias 
18S (long) (ingroup) | 1269 | 214(17%) |__ 374%) 0.54 1.16 towards transition substitutions. At 
, 18S (short) (ingroup) | 473 | 72.15%) | 802%) 0.79 146 | very extreme genetic distances, the 
28S (ingroup) | 583 |169(29%)| 40(7%) | 040 1.26 | transition-transversion ratio will 
Total (ingroup) | 1852 | 383 (21%) | _77(4%) | 0.34 1.19 | converge to 0.5. This is not apparent in 


percentage of parsimony-informative sites in both 
datasets and when outgroup taxa were removed, it 
also had the highest percentage of variable sites. 


The gamma shape parameter, which is an 
estimate of the rate variation across the sites in the 
alignment, ranged from 0.24-0.31 for alignments 
that included all taxa, and from 0.34-0.79 for 
ingroup alignments. The increased estimate of rate 
variation across sites (lower gamma value) in the 


18S (long) 


All taxa 


Topology 


Ingroup 


a Theonella sp1 
Theonella sp2 
Discodermia 
Corallistes 


Theonella sp1 
Theonella sp2 
Discodermia 


Corallistes 


Theonella sp1 
Theonella sp2 
Corallistes 
Discodermia 


All taxa 


any of the datasets in these analyses. 
The short 18S alignments generally 
produced trees with low amounts of resolution. 
Frequently none of the hypothesised phylogenies 
in Figure 3 were seen in the resulting bootstrap 
partition tables. The longer 18S alignments 
provided a greater amount of resolution, with 
bootstrap proportions sometimes becoming quite 
high. The 28S alignments were also quite well 
resolved and these also yielded high bootstrap 
values, particularly for the phylogeny that is 


285 


Ingroup 


Total 


Alltaxa__Ingroup 


18S (short) 


All taxa___Ingrou 


ioe) 
N 


mA] Ww = A Juju 
j ai li TP : 


© 
F ive) I 
© Cc © w © 
ive) al 


FIG. 3. Phylogenetic reconstructions of the relationships between Discodermia, Theonella and Corallistes. The 
first column (topology) indicates the topology that is under consideration and the open circle indicates the 
internal branch whose support is being assessed. The second column indicates the type of analysis that was 
undertaken: LD - LogDet; ML - Maximum Likelihood; P - Parsimony. The results are in four consecutive 
blocks according to the gene region that was used in the analysis. Within each block, thc left side indicates the 
bootstrap proportions for the alignment that included all taxa and the right side indicates the results when only 


the ingroup sequcnces were used. 


favoured by the morphological data (Fig. 3A). 
The combined sequence dataset alignments also 
displayed a reasonable amount of signal. 


DISCUSSION 


The analysis of sequence statistics showed that 
the 28S rRNA genc region provides the greatest 
amount of information per unit sequence length 
(Table 4). Despite the fact that the 28S alignment 
was hampered by the necessity of removing a 
large hypervariable portion, it still contained a 
high number and percentage of variable and 
parsimony-informative sites. The 18S gene 
behaved in a less efficient way. The combined 
dataset always contained the largest number of 
variable and parsimony-informative sites, but it 
did not contain the highest percentage of these 
sites in any of the analyses. 


It is curious to note the behaviour of the 
estimated shape of the gamma parameter, a. The 
gamma shape parameter is an cstimation of rate 
variation across sites. Lower gamma parameters 
indicate a more severe amount of rate variation 
whilst higher numbers indicate that the evolut- 
ionary ratc is more equivalcnt at all sites. There 
was an obvious difference in values between the 
taxon-inclusion sets. When all taxa were 
considered, the gamma shape puramcter was 
always lower than when only ingroup taxa were 
analysed. The addition of more taxa is simply 
increasing the amount of variability at sites that 
are already free to vary. Conserved sites remain 
unchanged with the addition of more taxa. 
Although the addition of more taxa has the effect 
of increasing the percentage of variable and 
parsimony-informative sites, the gamma shape of 
rate variation across sites is more marked. 


Not all of the three phylogenetic trees were 
congruent with the morphological phylogeny. 
The shorter of the two 18S rRNA alignments was 
unable to resolve the relationships of the four taxa 
of interest with any degree of confidence. Indeed 
during some bootstrap replicates, some other 
topologies, not considered in Figure 3 were found. 
The main reason for this was a complete lack of 
variability in the dataset. 


The longer 18S gene region was slightly more 
decisive about branching order. The topology 
that received strong support using this region was 
a pairing of Discodermia and Corallistes to the 
exclusion of the other taxa. However, this high 
level of support was only achieved using parsi- 
mony tree reconstruction. Given the general lack 
of confidence using the other methods, and the 
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inability of parsimony to compensate for super- 
imposed substitutions, it is possible that these 
findings are a result of long branch artifacts. This 
clade is rooted by a particularly long branch 
leading to the other demosponge taxa. The topol- 
ogy that was observed the least number of times 
during bootstrapping was the pairing of 
Discodermia with Theonella, and this was irresp- 
ective of type of analysis or taxon-inclusion set. 


The results for the 28S gene sequences were 
considerably different to those seen in the 18S 
analyscs. For this gene, the placement of 
Corallistes as the sister taxon to Theonella was 
never seen in any analysis (Fig, 3C) Of the two 
remaining alternative topologies, the placement 
of Discodermia as sister taxon to Theonella (Fig. 
3A) received considerably more support than the 
placement of Discodermia with Corallistes (Fig 
3B). Reasonably high bootstrap support was seen 
for Discodermia and Theonella as sister-groups 
using all methods of analyses irrespective of 
whcther the outgroup taxa were used or not. 


The alignment combining 18S and 28S rDNA 
data also yielded ambiguous results. The 
topology that suggests a sister taxon relationship 
between Corallistes and Theonella is very poorly 
supported (Fig. 3C). The other two topologies are 
more strongly supported, but neither was supp- 
orted with any degree of confidence and the 
differences in the levels of support do not justify 
acceptance of one over the other. It is likely that 
the ambiguous nature of the results from the 18S 
rRNA gene have a detrimental cffect on the 
combined alignment. During bootstrapping, a 
character can be selected from either gene region. 
Given that the 18S gene region is approximately 
150% larger than the 28S rRNA gene region, it 
will probably contribute more to each replicate 
on the whole. The result of this seems to be the 
carry-over of the ambiguous results from the 
separate analysis that used only the 18S rRNA 
gene. 


The topology that was most strongly supported 
using the 285 gene is consistent with the 
hypotheses of relationships deduced from 
morphological characters; Discodermia and 
Theonella are more closely related to each other 
than they are to Corallistes, and they are the more 
derived taxa. The 285 rRNA gene generally has a 
higher proportion of variable and 
parsimony-informative sites and can provide the 
best possibility of resolving porifcran phylo- 
genetic relationships, at least at the sub-ordinal 
level. 


RESOLUTION POTENTIAL OF 18S AND 28S rRNA SPONGE GENES 
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